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We show that in hydrodynamic simulations for relativistic heavy-ion coUisions, strong resonance 
decay calculations can be performed with fewer species of particle resonances while preserving good 
accuracy in single particle spectra and flow anisotropics. Such partial resonance calculations boost 
computation efficiency by a factor of 10 which is essential for large scale event-by-event simulations. 
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I. INTRODUCTION 



Over the past decade, relativistic hydrodynamics has 
established itself as an indispensable component in mod- 
eling the collective dynamics of the quark-gluon plasma 
(QGP) produced in relativistic heavy ion collisions 
l2l| . The properties of the QGP, in particular its 
shear viscosity, can be deduced by studying its effi- 
ciency of converting initial geometric anisotropics into 
final anisotropic flows. Because of this, measurements 
of collective flow anisotropics d„ (initially the elliptic 
flow fo, more recently higher-order harmonic flows w„>3 
|2a - l27t ) have generated strong interest and fertilized the 
study of initial geometry fluctuations and their dynami- 
cal consequences |28l - l4^ . 

The tool of choice for such studies is event-by-event 
hydrodynamics [MS, IHM EHU where each fluctu- 
ating initial density profile is evolved separately, followed 
by taking an average over the event ensemble to compute 
PT-spectra, anisotropic flow coefflcients and two-particle 
correlations in the flnal state which can then be compared 
with experimental data. Due to the limited number of fi- 
nal state particles in each event, these observables can be 
measured with good statistical precision only as ensem- 
ble averages. Due to algorithmic progress over the last 
few years, the hydrodynamic evolution part is no longer 
the bottleneck in such event-by-event studies; at least for 
(2-f-l)-dimensional simulations (which assume longitudi- 
nal boost-invariance) the largest fraction of the computer 
time is spent converting the hydrodynamic output into 
final particle distributions, cither on a "switching sur- 
face" between a macroscopic hydrodynamic description 
of the QGP fluid and a microscopic kinetic evolution of 
the dilute late hadronic rescattering stage 0, [s^l , or on a 
"kinetic decoupling" surface marking the transition from 
a strongly coupled fluid directly to a non-interacting gas 
of free-streaming hadrons. The high numerical cost of 
this "hydro-to-particle conversion" process results from 
the large number of unstable hadron resonances that need 
to be included and whose post-freeze-out decays (mostly 
due to strong-interaction processes, although for some 
comparisons with experimental data that have not been 
corrected [Sll, [IH for weak-decay feed-down, weak and 
electromagnetic decays must also be considered) modify 
the flnally observed particle distributions. 



The hydro-to-hadron conversion algorithm is based on 
the Cooper-Frye formula [53j which expresses the final 
hadron momentum distribution as an integral of the lo- 
cal equilibrium (for ideal fluid dynamics) or slightly off- 
equilibrium (in viscous fluids) distribution function for 
the particle species in question over the conversion sur- 
face. Contributions to the spectra of experimentally mea- 
sured stable particles from the strong decays of unstable 
resonances are then calculated from the single-particle 
spectra for the resonances [s^l- This requires the calcu- 
lation of the directly emitted ("thermal") particle mo- 
mentum distributions for all ^ 300 hadron species with 
mass typically up to 2 GeV via Cooper-Frye integrals, fol- 
lowed by the evaluation of the phase-space integrals [s^l 
for all contributing decay channels. On a typical personal 
computer with a single-core CPU this calculation takes 
presently about 2-3 hours, compared to 10-15 minutes 
for the preceding hydrodynamic evolution with, say, the 
(2-|-l)-dimensional viscous fluid code VISH2+1 [9l.[To|. 

The ^ 2 GeV cutoff in resonance mass is dictated by 
requiring convergence of the relative particle yields of 
the measured hadronic final state after all unstable res- 
onances have been allowed to decay. (The pion yields 
are especially sensitive to resonance feeddown.) Ex- 
perimental evidence points to chemical decoupling at 
a temperature of Tchom ~ 165 MeV, i.e. close to the 
(pseudo)critical temperature for the quark-hadron phase 
transition [ssj : at this temperature, only resonances with 
masses above 2 GeV are sufficiently strongly Boltzmann- 
suppressed that their decay contributions to stable par- 
ticle yields can be safely ignored. 

Here wc show that for an accurate determination of 
the pion and proton anisotropic flow coefficients w„ a 
much smaller number of resonances needs to be taken 
into account than for the hadron yields, and that even 
the shape of the azimuthally averaged pion and proton 
transverse momentum spectra can be reliably determined 
by accounting for only a small subset of the ~ 300 res- 
onance species mentioned above. These are the observ- 
ables needed for an extraction of the QGP shear viscosity 
from heavy- ion collision experiments . By rearranging 
the resonance decay table in order of decreasing impor- 
tance for the calculation of pT-spectra and coefficients 
instead of increasing mass, good convergence for these 
observables can be achieved with a significantly reduced 
set of only about 20-30 resonances. This speeds up the 
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computation by a factor of 10 - a significant gain in effi- 
ciency for the iterative determination of the QGP shear 
viscosity. 

The analysis presented here uses final states generated 
with the (2+l)-dimensional boost-invariant viscous hy- 
drodynamic code VISH2+1 for 200 ^GeV Au-|-Au colli- 
sions at the Rclativistic Heavy-Ion Collider (RHIC) and 
for 2.76 ^TeV Pb+Pb collisions at the Large Hadron 
Collider (LHC) at various collision centralities, with pre- 
viously determined [ill, E^, hydrodynamic input pa- 
rameters. We find very similar results at both collision 
energies and therefore show here only plots for LHC col- 
lisions. Since the decay contributions from different reso- 
nances to the mentioned observables depend only on their 
decay channels and transverse momentum distributions, 
we expect little sensitivity to the assumption of longitudi- 
nal boost-invariance implicit in our approach and expect 
our reordered resonance decay tables to perform equally 
well for both (2-|~l)-d and (3-|-l)-d hydrodynamic simula- 
tions, and for a wide range of input parameters (such as 
QGP viscosity, thermalization time, initial entropy and 
energy density, etc.). 



II. RESONANCE ORDERING 

The momentum distributions of directly emitted 
("thermal") resonances of species i are computed from 
the Cooper- Frye formula [Sj] : 



E 



dN 



9i 



d^p dy pt dpT d(j)p (2tt)'^ 



p''d^<Jf.iU + Sh). 



(1) 

Here E is the hydro-to-hadron conversion hypersurface, 
d^CT^ its surface normal vector, fio = l/[e'^(P "~''') q= 1] is 
the Bose or Fermi thermal equilibrium distribution func- 
tion, and Sfi accounts for viscous corrections (driven by 
the viscous pressure tensor tt^^ (x) on the conversion sur- 
face) of the local phase-space distribution along E. We 
assume the quadratic form Sf = i/o(l ± /o)- 



2J'J\^ jyi T2 (e+p) ■ 

Resonance decays increase the total yields of the sta- 
ble hadrons and change their momentum distributions. 
For kinematic reasons, most of the light decay daughters 
have low transverse momenta, thus modifying the shape 
of light stable hadrons (pions, kaons) particle spectra 
mostly in the region px < 1.5 GeV [5J]. We denote the 
total decay contribution to the momentum distribution 
of stable hadron species i by S[dNi/ (dyd'^pT)) , and the 
total spectrum (obtained by adding this to the thermally 
emitted spectrum dNf'^ / [dyd'^pr)) by dNl°^ / [dyd^pr)- 
(We here include only strong and electromagnetic de- 
cays.) The pT-integrated total yield S{dNi/dy) of de- 
cay products of species i is denoted by 6Ni, with 
7V*°* = m'^+SN, = N^^+ V, L^.m'^, where the sum is 



over resonances j and bj^i is the effective branching ratio 
(see Eq. (jH) below) for the decay j ^i. 

The contribution to 5Ni from a particular resonance 
J is not only influenced by its mass (through the Boltz- 



mann suppression factor e ^^^'^), but also by its spin 
degeneracy factor gj and its branching ratio bj^t into 
the decay channel that feeds stable particle species i. 
For each stable hadron species i it is therefore a dif- 
ferent set of resonances that makes the most important 
contributions. Our goal is to order the resonances in 
decreasing order of importance, for each stable particle 
species i. We here assume that the conversion surface 
has constant temperature Tconv The different hadron 
resonances have Tconv-dependent non-equilibrium fugac- 
ities Xj that ensure constant stable particle ratios equal 
to their chemical equilibrium values at Tchem and /is = 0, 
independent of the hydro-to-hadron conversion temper- 
ature Tconv While the actual fractions contributed by 
each resonance to the stable particle yields depend on 
Tconv, the ordering of these fractions is largely Tconv- 
independent. 

We start from the resonance table in the AZHYDRO pack- 
agcQ which includes 319 species of hadrons (counting 
different isospin states such as 7r+, tt*', Tr~ as separate 
species) with rest masses up to 2.25 GeV. After fixing 
the value of Tconv we look up the non-equilibrium fu- 
gacity Xj for each of these 319 species from the EOS 
s95p-PCE tables constructed in Ref. [5^. For each sta- 
ble particle species i we then generate an ordered list of 
resonances j that can decay directly into i. Note that in 
this ordering we account not only for direct decay con- 
tribution but also for multi-step decay cascades, where j 
first decays into an unstable resonance k which further 
decays (directly or through more intermediate steps) into 
the stable species i. 

Table |I] shows the beginning of this contribution table 
for positively charged pions, for a conversion temperature 
Tconv = 120 MeV. The "total contribution" percentages 
Cj^i in the third column are computed as 
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(3) 



where the effective branching ratios bj^i in Eq. ^ ac- 
count for multi-step decay cascades as follows: 



■ ^ bj^kibki^k2bk2 

ki,k2 



(4) 



The sum over in ([3]) takes care of quantum statistical 
effects, with the upper (lower) sign for bosons (fermions). 



1 AZHYDRO is available at the URL 

^http : //www . physics . ohio-state ■ edu/~f roderma/ | 
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FIG. 1. (Color online) Transverse momentum spectra of tt^, for a bumpy central (0-10% centrality, top panels) and a smooth 
peripheral (60-70% centrality, bottom panels) Pb-Pb collision at LHC energies. Panels (a) and (c) present the absolutely 
normalized spectra, while panels (b) and (d) show the normalized ratio r,r(pT) defined in Eq. ([5]). Different lines correspond to 
different cumulative resonance decay contributions between 0% ("thermal") and 100%. See text for discussion. 



name 


mass (GeV) 


total contribution (%) 


LJ 


0.78260 


15.398 


P" 


0.77580 


11.179 


P+ 


0.77580 


11.098 









TABLE I. Example of the vr'*' contribution table for 
Tconv = 120 MeV. 



For all hadrons except pions accurate results can be ob- 
tained by keeping only the first term fc= 1, i.e. by ig- 
noring quantum statistical effects. Even for pions, a few 
/c-terms suffice for good precision (in our calculations we 
truncate the series in (jS]) at fc = 10) . The complete or- 
dered resonance decay contribution tables for tt"*", K~^, 



p. A, and are given in the Appendix. Horizontal 
lines in the tables indicate where the cumulative reso- 
nance decay contributions c™* = X]j=i ^j^i exceed cer- 
tain threshold percentages (as indicated) of the total res- 
onance decay contribution to species i. 

In the following section we show the stable hadron px- 
spcctra and their anisotropic flow coefficients as a func- 
tion of these cumulative decay contribution percentages 
c™', in order to assess how many resonances from these 
ordered decay tables should be included for an accurate 
computation of these observables. 



III. RESULTS AND DISCUSSION 

Using the ordered tables described in Sec. lU and trun- 
cating the sum over resonance decay contributions at jcut 
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FIG. 2. (Color online) Same as Fig[l] but for protons. 



values corresponding to various different cumulative res- 
onance decay contribution tliresliolds c^"*, we performed 
calculations for 7r+, K'^ , and p. We tested individ- 
ual bumpy as well as (ellipticity-aligned and ensemble- 
averaged) smooth initial conditions at both RHIC and 
LHC energies, for a variety of collision centralities. Since 
the results were found all to be qualitatively similar, we 
show only a small selection, focussing on pious and pro- 
tons from one bumpy Pb-Pb event from the 0—10% cen- 
trality class and from the smooth averaged initial condi- 
tion corresponding to the 60—70% centrality class, both 
at LHC energy (v^ = 2.76 AGeV). 

Figure [T] shows the pion p^-spectra, for the bumpy cen- 
tral collision in the upper panels and the smooth periph- 
eral event in the lower panels. The left panels show the 
usual semilogarithmic plots of the absolutely normalized 
PT-distribution. As is well-known, the directly emitted 
("thermal" ) pious constitute only about 50-60% of all ob- 
served pions, the rest coming from resonance decays. The 



"thermal" spectrum also has the wrong shape: resonance 
decay pions predominantly contribute to the low-py part 
of the spectrum, making it steeper. However, this shape 
difference between the truncated and full resonance de- 
cay spectrum disappears almost completely already when 
including only the 9 strongest decay channels, account- 
ing for just 60% of the total pion yield from resonance 
decays. This is shown in the right panels of Fig. [1] where 
we plot the ratio 

rfJVf , 1 y-ic„t dNl''> 

nipr) = ' n (5) 

dypTdpT ^-^J — 1 dypTdpT 

for i = 7r as a function of px- {N^'''' is the contribution 
to particle species i from decays of particle species j (see 
Eq. ([1])); and jmax is the index of the last resonance in 
the ordered resonance decay table from SecHIl) The nu- 
merator includes only resonance decays up to jcut; but 
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FIG. 3. (Color online) The differential elliptic (u2, left panel) and triangular flow (us, right panel), for tt^ (upper panels) and 
p (lower panels), for one bumpy Pb-Pb event from the 0—10% centrality class at LHC energy. As in Fig. [1] lines of different 
styles and colors correspond to different cumulative resonance decay fractions. 



we renormalizc those decay contributions by the cumu- 
lative decay contribution c™' corresponding to the same 
jcut value, (c™* is easily calculated from Eqs. (|2l3p and 
directly obtained by summing the entries in the third 
column of the resonance decay table.) This renormaliza- 
tion corrects for the missing yield from the truncation of 
the decay table. The remaining effect (after missing yield 
renormalization) of the truncation on the shape of the px- 
spectrum is seen in panels (b) and (d) of Fig. [T] Whereas 
without any resonance decays the ratio rT^ipr) changes 
by almost a factor 2 between pt = and 2 GeV, this vari- 
ation is reduced to less than 5% already for c™* = 60%, 
for both bumpy and smooth initial conditions in both 
central and peripheral collisions. 

In Fig. [2] we show in the same way the proton spec- 
tra. Again the shape of the spectra can be accurately 
reproduced by taking into account a small fraction of all 
decay contributions (note the expanded vertical scale in 
Figs.[2]3,d): after renormalization to account for the miss- 
ing yield, just the 4 strongest of 75 decay channels (three 
charge states of the A(1232) resonance and one charge 
state of A(1600)), corresponding to 60% of the total res- 
onance decay yield for protons, reproduce the full proton 
spectrum with < 5% error between px = and 2 GeV. 

We conclude that, by accounting for the missing yield 



through appropriate renormalization, the correctly nor- 
malized total pion and proton spectra can be obtained, 
with shape errors < 5%, by including only the strongest 
decay channels accounting for the leading 60% of the to- 
tal resonance decay yields. A quick look at the tables in 
the Appendix shows that this will reduce the number of 
resonance decays (and thus computer time) by at least a 
factor 10. 

We now proceed to a discussion of the differential and 
px-intcgrated anisotropic flow coefficients u„, defined by 

J dypTdpTd(j)p I dyprdpT 

e"'/'" ^ [pT dpT dcl^, e»^^ , f , , / ^ . (6) 
J dypTdprdcpp / dy 

Here the spectrum dN / {dypTdpTd4>p) includes all contri- 
butions from the ordered resonance decay tabic for the 
considered stable species up to a certain threshold jcut, 
with the truncated resonance decay contribution renor- 
malized for the missing yield by a factor 1 /c™* as shown 
in the numerator of Eq. ([5]). In Figs. [3] and |4] below 
we specify the cumulative decay contribution percentage 
^cut indicate the truncation level corresponding to each 
curve. 



6 



0.03 
0.025 
0.02 



fo.oisl =— 

0.01 
0.005 




7r+, 0-10%, bumpy event 



-thermal 

-60% 

-80% 

-90% 

-95% 

-98% 

-100% 



0.035 



?5 



123456789 

harmonic order n 



0.1 



(c) 



7r+, 60-70%, smooth event 



0.03 
0.025 

0.02^ 
0.015 

0.01 
0.005 


0.1 



p, 0-10%, bmiipy event 



123456789 

harmonic order n 



p, 60-70%, smooth event 



0.08 



0.08 



S 



0.06 
0.04 
0.02 



0.06 
0.04 
0.02 



123456789 123456789 

harmonic order n harmonic order n 

FIG. 4. (Color online) The pr-integrated anisotropic flow coeflicients Vn, n = l, . . . ,9, for tt^ (a,c) and p (b,d), for a bumpy 
central (0-10%) collision event (a,b) and a smooth averaged peripheral (60-70%) collision event (c,d). Line colors and styles as 
in Figs. IH3I For a discussion of the two sets of lines for each harmonic order n see text. 



Figure |3] shows the difFcrcntial elliptic and triangular 
Hows for pions and protons, for one single bumpy central 
(0-10% centrality) event. We see that once again excel- 
lent agreement with the full resonance decay calculation 
is already obtained when including only the small subset 
of resonances that account for the top 60% of the reso- 
nance decay yields. We checked that this result is generic, 
i.e. it does not depend on the selected event (although 
the elliptic and triangular flows do). 

For the pr-integrated harmonic flow coefficients w„ we 
show in Fig. 2] results for all harmonic orders from n = l 
to 9, again for pions and protons and for a bumpy cen- 
tral as well as a smooth peripheral event. For the smooth 
averaged initial condition, the odd harmonics vanish by 
symmetry. For fluctuating initial conditions, the v„ val- 
ues shown here and their relative size depend on the ran- 
domly selected event. All plots shown in this paper are 
based on one and the same bumpy central collision event. 

For each harmonic order n, Fig. |4] shows two sets of 
results. The left set corresponds to results obtained by 
using the truncated resonance decay spectra shown in 
Figs. [T^,c and [2^,c, without missing yield renormaliza- 
tion. The right set uses the renormalized truncated decay 
spectra as defined in the numerator of Eq. ([5]) . One ob- 
serves a much faster convergence towards the full result 



in the right sets than in the left sets. The reason is that, 
by renormalizing the truncated resonance decay contri- 
butions for the missing yield, the correct mixing ratio be- 
tween direct thermal and indirect decay contributions is 
ensured and the shape of the total py-spectrum is approx- 
imated much more accurately than without renormaliza- 
tion (see Figs. [Tl3,d and[5]3,d). Figure 0] demonstrates 
that, when using the renormalized truncated decay spec- 
tra, accounting for just the top 60% decay contributions 
(i.e. including only the 9 strongest decay channels con- 
tributing to pions and the 4 strongest decay channels 
contributing to the proton spectra) reproduces the full 
results for the harmonic flow coefficients w„ with excel- 
lent precision: The lines corresponding to different c™* 
values >60% are almost indistinguishable. 

Future precision extractions of the QGP viscosity may 
require highly precise u„ values. For such a purpose one 
can adjust jcut to include a larger fraction of all resonance 
decays if needed. 

For a given precision, the required minimal jcut trun- 
cation indices and cumulative resonance decay fractions 
c^' for kaons lie between those for pions and protons. 
The c^"* for i=Tr,K,p are almost identical at RHIC and 
LHC energies, i.e. only weakly sensitive to radial flow. 
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IV. CONCLUSIONS 

We have shown that for a sufficiently accurate deter- 
mination of the differential anisotropic flow coefficients 
VniPr) only those resonances need to be included that 
generate the top 60% of the largest decay contributions 
to the stable particle yields. For the single particle spec- 
tra, correct normalization of the total yield requires a 
renormalization of the truncated resonance decay yield 
as given in the numerator of Eq. ([5]). With this renormal- 
ization, good convergence of the slope of the pion spectra 
and of the pT-integrated anisotropic flow coefficients Vn 
requires inclusion of only the 9 strongest contributing 
channels for pious and only the 4 strongest channels for 
protons, accounting in both cases for just 60% of the total 
decay yield. This reduces the number of resonance decay 
channels to be evaluated by a factor >fO, without loss 
of precision, leading to a similar reduction of the total 



computing time for the final stable hadron distributions. 

In hybrid model calculations [s^l the late hadronic 
stage is described microscopically by a Boltzmann cas- 
cade that propagates a reduced set of resonances un- 
til final kinetic decoupling. In this case the spectra of 
all unstable resonances that are explicitly included in 
the Boltzmann cascade must be generated on the con- 
version surface. This is still only a small subset of all 
resonances included in the resonance decay tables. The 
optimal ordering of the resonance decay tables for the 
purpose of generating input for the late-stage Boltzmann 
cascade and the corresponding optimized truncation frac- 
tions c™' will be studied in a follow-up report. 
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V. APPENDIX: FEED DOWN CONTRIBUTION 
TABLES FOR %, K, p, A, E+, AND H 



name 


mass (GeV) 


total contribution (%) 






0.78259 


15.398 






0.7758 


11.179 




P+ 


0.7758 


11.098 




K*+{892) 


0.89166 


5.54 




K*°{892) 


0.8961 


5.355 




V 


0.54775 


4.682 




A"(1232) 


1.232 


2.613 




A++(1232) 


1.232 


2.606 




6+(1235) 


1.2295 


2.498 


60% 


r;'(958) 


0.95778 


2.069 




aj(980) 


0.9847 


1.862 




/ii(1170) 


1.17 


1.26 




(1260) 


1.23 


1.226 




6;- (1235) 


1.2295 


1.19 




b'i'(1235) 


1.2295 


1.181 




a+(1320) 


1.3183 


1.177 




E+(1385) 


1.3828 


1.09 




E''(1385) 


1.3872 


1.057 




/i(1285) 


1.2818 


0.994 




iff (1270) 


1.273 


0.963 




A'i*"(1270) 


1.273 


0.962 




A°(1232) 


1.232 


0.857 




A+(1232) 


1.232 


0.856 




a;(1260) 


1.23 


0.817 




a^(1320) 


1.3183 


0.771 




(^(1020) 


1.0195 


0.762 


80% 


/o(980) 


0.9741 


0.613 




iff (1270) 


1.273 


0.472 




Jf-(1270) 


1.273 


0.472 




/2(1270) 


1.2754 


0.45 




a" (1260) 


1.23 


0.409 




(980) 


0.9847 


0.402 




ag(980) 


0.9847 


0.399 




a^(1320) 


1.3183 


0.398 




7r+(1400) 


1.376 


0.373 




K*+{U30) 


1.4256 


0.368 




E°(1530) 


1.5318 


0.357 




A-(1600) 


1.6 


0.356 




A++(1600) 


1.6 


0.356 




7?|°(1430) 


1.4324 


0.35 




=-(1530) 


1.535 


0.347 




7?? (1400) 


1.402 


0.341 




K+iUOO) 


1.402 


0.34 




p(1520) 


1.52 


0.305 




n(1520) 


1.52 


0.305 




?7(1295) 


1.294 


0.297 




i?*-(1410) 


1.414 


0.294 




Jf*+(1410) 


1.414 


0.294 





TABLE II. Resonance contribution list for tt^, 
T.onv = 120 MeV. 



name 


mass (GeV) 


total contribution (%) 




7r+(1300) 


1.3 


0.257 




A''(1600) 


1.6 


0.249 




A+(1600) 


1.6 


0.249 




ri(1440) 


1.44 


0.241 




p(1440) 


1.44 


0.241 




a;(1420) 


1.419 


0.205 




aj(1450) 


1.474 


0.196 




A(1405) 


1.4065 


0.188 


90% 


A(1405) 


1.4065 


0.188 




J72(1645) 


1.617 


0.185 




/i(1420) 


1.4263 


0.174 




JC?(1400) 


1.402 


0.173 




K-{UOO) 


1.402 


0.173 




ag(1450) 


1.474 


0.153 




n(1675) 


1.675 


0.147 




p(1675) 


1.675 


0.147 




A-(1700) 


1.7 


0.145 




A++(1700) 


1.7 


0.145 




a;3(1670) 


1.667 


0.144 




7r°(1300) 


1.3 


0.142 




JC*"(1410) 


1.414 


0.14 




if*- (1410) 


1.414 


0.14 




A°(1600) 


1.6 


0.135 




A+(1600) 


1.6 


0.134 




p+(1690) 


1.6888 


0.123 




p(1680) 


1.685 


0.119 




n(1680) 


1.685 


0.119 




A°(1700) 


1.7 


0.117 




A+(1700) 


1.7 


0.117 




k;-{u30) 


1.4256 


0.115 




A- (1620) 


1.63 


0.113 




A++(1620) 


1.63 


0.113 




A(1520) 


1.5195 


0.112 




A(1520) 


1.5195 


0.112 




7^1° (1430) 


1.4324 


0.109 




?7(1405) 


1.4103 


0.104 




n(1535) 


1.535 


0.103 




p(1535) 


1.535 


0.102 




n(1700) 


1.7 


0.102 




p(1700) 


1.7 


0.102 




n(1720) 


1.72 


0.098 




p(1720) 


1.72 


0.098 




P§(1690) 


1.6888 


0.097 




7r-(1300) 


1.3 


0.093 




ao (1450) 


1.474 


0.091 




p(1520) 


1.52 


0.088 


95% 


ri(1520) 


1.52 


0.088 




A''(1620) 


1.63 


0.086 




A+(1620) 


1.63 


0.086 




7r2+(1670) 


1.6724 


0.086 





TABLE III. Resonance contribution list for 7r+, 
Tconv = 120 MeV (continued). 



name 


mass (GeV) 


total contribution (%) 


p°(1450) 


1.465 


0.082 


p+(1450) 


1.465 


0.082 


7rj7(1400) 


1.376 


0.08 


7r?(1400) 


1.376 


0.08 


E"(1670) 


1.67 


0.079 


E+(1670) 


1.67 


0.079 


E+(1775) 


1.775 


0.075 


I)"(1775) 


1.775 


0.075 


/9+(1700) 


1.72 


0.072 


n(1700) 


1.7 


0.072 


p(1700) 


1.7 


0.072 


n(1440) 


1.44 


0.072 


p(1440) 


1.44 


0.072 


A°(1700) 


1.7 


0.07 


A+(1700) 


1.7 


0.07 


A(1690) 


1.69 


0.069 


A(1690) 


1.69 


0.069 


S°(1385) 


1.3837 


0.069 


E°(1385) 


1.3837 


0.069 


cj(1650) 


1.67 


0.067 


n(1675) 


1.675 


0.061 


p(1675) 


1.675 


0.061 


7rj^(1600) 


1.653 


0.06 


^§(1770) 


1.773 


0.056 


2. V / 


1.773 


0.056 


E°(1670) 


1.67 


0.05 


E°(1670) 


1.67 


0.05 


A"(1905) 


1.89 


0.05 


A++(1905) 


1.89 


0.05 


^•^(1700) 


1.72 


0.049 


A''(1620) 


1.63 


0.048 


A+(1620) 


1.63 


0.048 


71-^(1670) 


1.6724 


0.045 


n(1710) 


1.71 


0.044 


A''(1905) 


1.89 


0.044 




1.412 


0.044 


A+(1905) 


1.89 


0.044 


p(1710) 


1.71 


0.044 


i?o''(1430) 


1.412 


0.044 


E"(1660) 


1.66 


0.043 


E+(1660) 


1.66 


0.043 


n(1650) 


1.655 


0.042 


p(1650) 


1.655 


0.042 


fo (15001 

J u V / 


1.507 


0.041 


7)(1475) 


1.476 


0.04 


A-(1950) 


1.93 


0.04 


A++(1950) 


1.93 


0.04 


A-(1920) 


1.92 


0.039 



TABLE IV. Resonance contribution list for tt 
Tconv = 120 MeV (continued). 



name 


mass (GeV) 


total contribution (%) 


A++(1920) 


1.92 


0.039 


A(1830) 


1.83 


0.037 


A(1830) 


1.83 


0.037 


p7(1690) 


1.6888 


0.037 


i?;°fl780) 


1.776 


0.034 


K*+(1780] 


1.776 


0.034 


p" (1700) 


1.72 


0.034 


n(1535) 


1.535 


0.033 


p(1535) 


1.535 


0.033 


n(1720) 


1.72 


0.033 




1.72 


0.033 


p(1680) 


1.685 


0.032 


n(1680) 


1.685 


0.032 


if? (1770') 


1.773 


0.032 


7?? (1820) 


1.816 


0.032 


7<',"(1770) 


1.773 


0.032 


Kt (1820) 


1.816 


0.032 


7r^(1670) 


1.6724 


0.031 


A(1600) 


1.6 


0.03 


A(1600) 


1.6 


0.03 


fo(1370) 


1.4 


0.029 


A''(1905) 


1.89 


0.028 


A+(1905) 


1.89 


0.028 


n(1710) 


1.71 


0.027 


p(1710) 


1.71 


0.027 


0(1680) 


1.68 


0.026 


A°(1950) 


1.93 


0.024 


A+(1950) 


1.93 


0.024 


ii'*+(1680) 


1.717 


0.024 


A'*'' (1680) 


1.717 


0.024 


7rj"(1600) 


1.653 


0.024 


7r?(1600) 


1.653 


0.024 


E"(1915) 


1.915 


0.023 


E+(1915) 


1.915 


0.023 


A°(1920) 


1.92 


0.023 


A+(1920) 


1.92 


0.023 


H°(1820) 


1.823 


0.023 


H"(1820) 


1.823 


0.023 


E~(1940) 


1.94 


0.022 


S+(1940) 


1.94 


0.022 


A(1670) 


1.67 


0.021 


A(1670) 


1.67 


0.021 


A"(1910) 


1.91 


0.02 


A++(1910) 


1.91 


0.02 


E°(1775) 


1.775 


0.02 


S°(1775) 


1.775 


0.02 


^2^(1820) 


1.816 


0.019 


A2~(1820) 


1.816 


0.019 


A-(1930) 


1.96 


0.017 


A++(1930) 


1.96 


0.017 



TABLE V. Resonance contribution list for tt 
Tconv = 120 MeV (continued). 



name 


mass (GeV) 


total contribution (%) 


/2(1525) 


1.525 


0.015 


S°(1660) 


1.66 


0.014 


E°(1660) 


1.66 


0.014 


S"(1750) 


1.75 


0.014 


E+(1750) 


1.75 


0.014 


E°(1750) 


1.75 


0.014 


E°(1750) 


1.75 


0.014 


is:3*°(1780) 


1.776 


0.014 


is-*- (1780) 


1.776 


0.014 


/2(2010) 


2.011 


0.012 


A" (1600) 


1.6 


0.012 


A++(1600) 


1.6 


0.012 


A(1890) 


1.89 


0.011 


A(1890) 


1.89 


0.011 


A"(1920) 


1.92 


0.011 


A+(1920) 


1.92 


0.011 


A"(1910) 


1.91 


0.011 


A+(1910) 


1.91 


0.011 


A''(1950) 


1.93 


0.011 


A+(1950) 


1.93 


0.011 


A(1820) 


1.82 


0.011 


A(1820) 


1.82 


0.011 


7r+(1800) 


1.812 


0.01 


A(1800) 


1.8 


0.01 


A(1800) 


1.8 


0.01 


E°(1940) 


1.94 


0.009 


E°(1940) 


1.94 


0.009 


E-(1750) 


1.75 


0.009 


E+(1750) 


1.75 


0.009 


E°(1915) 


1.915 


0.009 


E0(1915) 


1.915 


0.009 


H°(1820) 


1.823 


0.009 


=-(1820) 


1.823 


0.009 


E°(1950) 


1.95 


0.008 


=-(1950) 


1.95 


0.008 


A(1810) 


1.81 


0.008 


A(1810) 


1.81 


0.008 


if*- (1680) 


1.717 


0.007 


A'*°(1680) 


1.717 


0.007 


E+(1775) 


1.775 


0.007 


E-(1775) 


1.775 


0.006 


A°(1930) 


1.96 


0.006 


A+(1930) 


1.96 


0.006 


r2(2250) 


2.252 


0.006 


A°(1910) 


1.91 


0.005 


A+(1910) 


1.91 


0.005 


n(1650) 


1.655 


0.005 


p(1650) 


1.655 


0.005 



TABLE VI. Resonance contribution list for n 
Tconv = 120 MeV (continued). 
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name 


mass (GeV) 


total contribution (%) 


7r"(1800) 


1.812 


0.004 


7r°(1800) 


1.812 


0.004 


(^3(1850) 


1.854 


0.004 


f2(1950) 


1.945 


0.003 


A" (1920) 


1.92 


0.003 


A++(1920) 


1.92 


0.003 


fofl710) 


1.715 


0.002 


A~(1910) 


1.91 


0.002 


A++(1910) 


1.91 


0.002 


E-(1940) 


1.94 


0.001 


£+(1940) 


1.94 


0.001 


E" (1915) 


1.915 





E+(1915) 


1.915 






TABLE VII. Resonance contribution list for 7r+, for 
Tconv = 120 MeV (continued). 



name 


mass (GeV) 


total contribution (%) 




K*°{892) 


0.8961 


35.857 




7!"+ (892) 


0.89166 


18.52 




(^(1020) 


1.0195 


16.036 


60% 


if+(1270) 


1.273 


3.631 




Kl (1270) 


1.273 


3.287 




a+(980) 


0.9847 


1.807 




/i(1420) 


1.4263 


1.754 


80% 


7^2*" (1430) 


1.4324 


1.61 




K+{UOO) 


1.402 


1.446 




K*+iim 


1.4256 


1.241 






1.414 


1.217 




is'? (1400) 


1.402 


1.127 




if*°(1410) 


1.414 


1.033 




a!^(980) 


0.9847 


0.898 




/o(980) 


0.9741 


0.869 


90% 


A(1520) 


1.5195 


0.798 




/2(1525) 


1.525 


0.711 




/i(1285) 


1.2818 


0.493 




=-(1690) 


1.69 


0.462 




E-(1775) 


1.775 


0.431 




7?(1475) 


1.476 


0.405 




A'o*°(1430) 


1.412 


0.295 




E-(1820) 


1.823 


0.288 




aj(1320) 


1.3183 


0.274 




?7(1405) 


1.4103 


0.265 




(^(1680) 


1.68 


0.265 


95% 


/2(2010) 


2.011 


0.259 




E°(1775) 


1.775 


0.237 




E°(1690) 


1.69 


0.23 




A(1820) 


1.82 


0.196 




ii'3+(1780) 


1.776 


0.191 




^2+ (1770) 


1.773 


0.189 




E-(1670) 


1.67 


0.177 




£-(1750) 


1.75 


0.171 




iS'2°(1770) 


1.773 


0.166 




£-(1660) 


1.66 


0.165 




A(1600) 


1.6 


0.161 




?7(1295) 


1.294 


0.157 




ifo*+(1430) 


1.412 


0.148 




if^(1820) 


1.816 


0.145 




A(1690) 


1.69 


0.144 




A'3*°(1780) 


1.776 


0.14 




a§(1320) 


1.3183 


0.137 


98% 


if*" (1680) 


1.717 


0.117 





TABLE VIII. Resonance contribution list for K 
Tconv = 120 MeV. 



name 


mass (GeV) 


total contribution (%) 


(/)3(1850) 


1.854 


0.114 


f2(1270) 


1.2754 


0.099 


A(1810) 


1.81 


0.091 


E°(1670) 


1.67 


0.089 


A(1890) 


1.89 


0.087 


E°(1750) 


1.75 


0.085 


E°(1660) 


1.66 


0.082 




1.816 


0.082 


=-(1950) 


1.95 


0.082 


A(1670) 


1.67 


0.082 


i<r*+(1680) 


1.717 


0.078 


A(1800) 


1.8 


0.076 


p(1720) 


1.72 


0.072 


p(1710) 


1.71 


0.068 


??2(1645) 


1.617 


0.06 


H''(1820) 


1.823 


0.052 


a+(1450) 


1.474 


0.052 


£+(1940) 


1.94 


0.047 


E°(1940) 


1.94 


0.045 


n(2250) 


2.252 


0.044 


E+(1775) 


1.775 


0.043 


E-(1940) 


1.94 


0.043 


E-(1915) 


1.915 


0.036 


p(1650) 


1.655 


0.032 


ao(1450) 


1.474 


0.026 


fo(1710) 


1.715 


0.024 


/9t(1690) 


1.6888 


0.023 


E°(1915) 


1.915 


0.018 


A(1830) 


1.83 


0.017 


n?(1690) 


1.6888 


0.016 


7rJ(1670) 


1.6724 


0.015 


7r9d670~) 


1.6724 


0.013 


fo(1500) 


1.507 


0.012 


^^(1670) 


1.6724 


0.01 


fnfl370~) 


1.4 


0.01 


07(1690) 


1.6888 


0.008 


A++(1920) 


1.92 


0.004 


A+(1920) 


1.92 


0.003 


A++(1950) 


1.93 


0.003 


A+(1950) 


1.93 


0.002 


is:? (1400) 


1.402 


0.002 


i!r]"(1400) 


1.402 


0.002 


A°(1920) 


1.92 


0.001 


/2(1950) 


1.945 


0.001 


A°(1950) 


1.93 


0.001 


i?§(1820) 


1.816 





/sTa" (1820) 


1.816 






TABLE IX. Resonance contribution list for 
Tconv = 120 MeV (continued). 



A'+, for 



name 


mass (GeV) 


total contribution (%) 




A++(1232) 


1.232 


29.842 




A+(1232) 


1.232 


19.816 




A''(1232) 


1.232 


9.813 




A++(1600) 


1.6 


2.787 


60% 


n(1520) 


1.52 


2.487 




p(1520) 


1.52 


2.169 




A+(1600) 


1.6 


2.049 




p(1440) 


1.44 


2.034 




n(1440) 


1.44 


1.943 




p(1535) 


1.535 


1.452 




A++(1700) 


1.7 


1.386 




A(1520) 


1.5195 


1.365 




p(1675) 


1.675 


1.347 




A°(1600) 


1.6 


1.314 




p(1700) 


1.7 


1.228 


80% 


A++(1620) 


1.63 


1.135 




n(1680) 


1.685 


1.003 




A+(1700) 


1.7 


0.987 




n(1675) 


1.675 


0.975 




p(1680) 


1.685 


0.917 




n(1535) 


1.535 


0.793 




A+(1620) 


1.63 


0.785 




n(1720) 


1.72 


0.752 




E+(1775) 


1.775 


0.738 




A" (1600) 


1.6 


0.574 




A°(1700) 


1.7 


0.54 


90% 


A++(1905) 


1.89 


0.536 




p(1720) 


1.72 


0.487 




p(1710) 


1.71 


0.456 




A''(1620) 


1.63 


0.433 




n(1650) 


1.655 


0.429 




A++(1950) 


1.93 


0.411 




E°(1775) 


1.775 


0.405 




A+(1905) 


1.89 


0.366 




n(1700) 


1.7 


0.342 




A(1820) 


1.82 


0.334 




E+(1670) 


1.67 


0.302 




p(1650) 


1.655 


0.29 


95% 


E+(1750) 


1.75 


0.289 




E+(1660) 


1.66 


0.282 




A+(1950) 


1.93 


0.281 




A(1600) 


1.6 


0.275 




A(1690) 


1.69 


0.247 




A++(1920) 


1.92 


0.227 




n(1710) 


1.71 


0.203 




A++(1930) 


1.96 


0.196 




A°(1905) 


1.89 


0.189 




A+(1920) 


1.92 


0.182 





TABLE X. Resonance contribution list for p, 
Tconv = 120 MeV. 



name 


mass (GeV) 


total contribution (%) 




A(1810) 


1.81 


0.156 




A" (1950) 


1.93 


0.153 




S°(1670) 


1.67 


0.151 




A(1890) 


1.89 


0.15 


98% 


E''(1750) 


1.75 


0.145 




E°(1660) 


1.66 


0.141 




A(1670) 


1.67 


0.139 




S+(1940) 


1.94 


0.138 




A''(1920) 


1.92 


0.137 




A" (1700) 


1.7 


0.136 




A+(1930) 


1.96 


0.13 




A(1800) 


1.8 


0.13 


99% 


A++(1910) 


1.91 


0.121 




A+(1910) 


1.91 


0.096 




A-(1920) 


1.92 


0.092 




A"(1620) 


1.63 


0.081 




E°{1M0) 


1.94 


0.077 




E-(1775) 


1.775 


0.074 




A°(1910) 


1.91 


0.071 




A°(1930) 


1.96 


0.065 




E+(1915) 


1.915 


0.062 




A-(1910) 


1.91 


0.046 




S°(1915) 


1.915 


0.031 




A(1830) 


1.83 


0.029 




A- (1950) 


1.93 


0.023 




A- (1905) 


1.89 


0.016 




E-(1940) 


1.94 


0.015 


100% 



TABLE XI. Resonance contribution list for p, for 
Tconv = 120 MeV (continued). 
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name 


mass (GeV) 


total contribution (%) 




E« 


1.1926 


24.775 




E+(1385) 


1.3828 


17.893 




E-(1385) 


1.3872 


17.346 


60% 


E°(1385) 


1.3837 


16.555 




A(1405) 


1.4065 


3.103 




A(1520) 


1.5195 


2.18 


80% 


E-(1670) 


1.67 


1.304 




S+(1670) 


1.67 


1.297 




A(1690) 


1.69 


1.178 




H~(1690) 


1.69 


1.132 




H»(1690) 


1.69 


1.128 




E+(1775) 


1.775 


0.771 




E~(1775) 


1.775 


0.769 




A(1830) 


1.83 


0.754 


90% 


E°(1775) 


1.775 


0.749 




E-(1660) 


1.66 


0.708 




E+(1660) 


1.66 


0.707 




H°(1820) 


1.823 


0.707 




E:-(1820) 


1.823 


0.705 




E°(1750) 


1.75 


0.619 




A(1670) 


1.67 


0.607 




A(1600) 


1.6 


0.486 


95% 


E"(1670) 


1.67 


0.485 




E''(1660) 


1.66 


0.471 




E-(1915) 


1.915 


0.382 




E+(1915) 


1.915 


0.381 




E°(1915) 


1.915 


0.264 




A(1820) 


1.82 


0.216 




H°(1950) 


1.95 


0.202 




e:-(1950) 


1.95 


0.201 


98% 


E-(1940) 


1.94 


0.189 




E+(1940) 


1.94 


0.188 




n(1720) 


1.72 


0.176 




p(1720) 


1.72 


0.176 




A(1890) 


1.89 


0.174 




n(1710) 


1.71 


0.167 


99% 


p(1710) 


1.71 


0.167 




E°(1940) 


1.94 


0.155 




A(1800) 


1.8 


0.155 




n(1650) 


1.655 


0.078 




p(1650) 


1.655 


0.078 




E-(1750) 


1.75 


0.076 




E+(1750) 


1.75 


0.075 




A(1810) 


1.81 


0.05 




A''(1920) 


1.92 


0.007 




A+(1920) 


1.92 


0.007 




A°(1950) 


1.93 


0.004 




A+(1950) 


1.93 


0.004 


100% 



TABLE XII. Resonance contribution list for A, for 
Tconv = 120 MeV. 
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name 


mass (GeV) 


total contribution (%) 




A(1405) 


1.4065 


25.159 




A(1520) 


1.5195 


10.121 




E+(1385) 


1.3828 


9.288 




E°(1385) 


1.3837 


9.179 




E°(1670) 


1.67 


6.648 


60% 


E+(1670) 


1.67 


6.618 




H°(1690) 


1.69 


4.578 




E+(1750) 


1.75 


4.425 




A(1600) 


1.6 


3.969 




A(1690) 


1.69 


3.845 


80% 


A(1830) 


1.83 


2.243 




A(1670) 


1.67 


1.926 




E+(1660) 


1.66 


1.916 




E°(1660) 


1.66 


1.916 


90% 


H°(1820) 


1.823 


1.045 




E+(1915) 


1.915 


0.974 




E''(1915) 


1.915 


0.968 




E+(1775) 


1.775 


0.859 


95% 


E''(1775) 


1.775 


0.805 




A(1820) 


1.82 


0.791 




E~(1775) 


1.775 


0.6 




E+(1940) 


1.94 


0.405 


98% 


A(1810) 


1.81 


0.402 




E°(1940) 


1.94 


0.387 


99% 


A(1800) 


1.8 


0.377 




A(1890) 


1.89 


0.231 




E-(1940) 


1.94 


0.101 




A++(1920) 


1.92 


0.08 




A++(1950) 


1.93 


0.051 




A+(1920) 


1.92 


0.027 




S-(1750) 


1.75 


0.02 




E°(1750) 


1.75 


0.02 




A+(1950) 


1.93 


0.015 




E-(1915) 


1.915 


0.007 


100% 



TABLE XIII. Resonance contribution list for E^, for 
Tconv = 120 MeV. 



name 


mass (GeV) 


total contribution (%) 




S"(1530) 


1.5318 


62.049 


60% 


=-(1530) 


1.535 


30.103 


90% 


2°(1820) 


1.823 


2.42 




=-(1820) 


1.823 


2.334 


95% 


H"(1950) 


1.95 


1.427 


98% 


^(2250) 


2.252 


0.957 


99% 


=-(1950) 


1.95 


0.71 


100% 



TABLE XIV. Resonance contribution list for H , for 
Tconv = 120 MeV. 



